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Abstract 

We investigate the effect of friction on domain formation in disoriented chi- 
ral condensate. Including a friction term, we solve the equation of motion 
of the linear sigma model fields, in the Hartree approximation. With boost- 
invariance and cylinderical symmetry, irrespective of friction, on average, we 
donot find any indication of domain like formation with quenched initial con- 
dition. However, with or without friction, some events can be found with 
large instabilities, indicating possible DCC domain formation in those events. 
With friction time scale during which instabilities grows increases. Corre- 
spondingly, with friction, it is possible to obtain large sized domains in some 
particular events. 
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The possibility of forming disoriented chiral condensate (DCC) in relativistic heavy ion 
colhsions has generated considerable research activities in recent years. The idea was first 
proposed by Rajagopal and Wilczek They argued that for a second order chiral phase 

transition, the chiral condensate can become temporarily disoriented in the nonequilibrium 
conditions encountered in heavy ion collisions. As the temperature drops below Tc, the chiral 
symmetry begins to break by developing domains in which the chiral field is misaligned 
from its true vacuum value. The misaligned condensate has the same quark content and 
quantum numbers as do pions and essentially constitute a classical pion field. The system 
will finally relaxes to the true vacuum and in the process can emit coherent pions. Since the 
disoriented domains have well defined isospin orientation, the associated pions can exhibit 
novel centauro-like fluctuations of neutral and charged pions . 



Most dynamical studies of DCC have been based on the linear sigma model, in which 
the chiral degrees of freedom are described by the real 0(4) field $ = (c, 11), having the 
equation of motion, 

[n + X{^^ ~ v^)]^ = Hn, (1) 

The parameters of the model can be fixed by specifying the pion decay constant, /7r=92 
MeV and the meson masses, 771,^=135 MeV and mo-=600 MeV, leading to A = (m^ — 
m2)/2/2=20.14 and v = [{ml - 3ml)/ {ml -mD^/^f ^=86.71 MeV and H = (120.55Me\/)3 
|TB| . It is apparent from eq.|T] that the vacuum is aligned in the a direction ^^ac = {fw, 0) 
and at low temperature the fluctuations represent nearly free a and tt mesons. At very high 
temperature well above v, the field fluctuations are centered near zero and approximate 0(4) 
symmetry prevails. 

It is instructive to decompose the chiral field, 

$(r, t) =< 0(r, t) > +6(p{r, t) (2) 

where < >is the mean field and 5(f) are the semiclassical fluctuations around < (j) > and 
can be identified with quasi-particle excitations. Using eq.|^ and taking the average of eq.|l|, 
the equation of motion for the mean fields in the Hartree approximation can be obtained as 



— -V^ <(j)>= X{v^- < (f) >^ -3 < 6(j)l >-< 6(j)l >)<(/) > +Hn, (3) 

where < (p >^=< 0j >< 0j >, is the component of the fluctuation parallel to < > 
and 5(j)± is the orthogonal component. This equation imply that the motion of the mean 
field is determined by the effective potential, 

V{< (j)>) = ^{< (j) +3 < 6(j)l > + < S(f)l > -v"") (4) 

which clearly differs from the zero temperature one in presence of fluctuations. By varying 
the fluctuations, chiral symmetry can be restored or spontaneously broken. It is also evident 
that the evolution of the mean field critically depends on the initial values of the fluctuations. 
When 5^ = (3 < 50| > — < 6(j)']_ >)/6 is large enough the chiral symmetry is approximately 
(as Ht^ 0) restored. If the explicit chiral symmetry breaking term is neglected, the phase 
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transition takes place at the critical fluctuations 51 = For < S^, the effective 

potential takes its minimum value at < >= ((Te,0), where depends on 5^. When the 
mean flelds are displaced from this equilibrium point to the central lump of the Mexican hat 
(< >~ 0), the effective mass square 



' - X{v^- < (f) >^ -3 < 6(j)l > - < S(j)l >) (5) 



will become negative and DCC can form. Since the domain size is directly related to the 
time scale, during which the effective mass remains negative, it strongly depends on the 
initial condition of the system. By varying the < (pi > and 6^, quench or annelaing like 
initial condition can be obtained . In an important paper Asakawa et al studied eq.|^ 



with initial condition corresponding to quench and annealing. They found that domains of 
disoriented chiral condensate with 4-5 fm in size can from through a quench. Annealing on 
the otherhand, leads to smaller sized domains. 

In the present paper, we solve eq.|^ with a diffrent motivation. Our interest is to investi- 
gate the effect of friction on DCC domain formation. Dissipative effects like friction damp 
the motion of the flelds, inhibiting large oscillations. Naively, one would expect reduction of 
DCC domain formation with friction. Recently, Biro and Greiner using the Langevin 



equation for the linear sigma model, investigated the interplay of friction and white noise on 
the evolution of the order parameter. While noise greatly diminishes the possibilty of DCC 
domain formation, in some orbits, large instabilities can result, producing DCC domains. 
We have also studied the effect of friction on DCC domain formation using the Langevin 



equation [|I6|. There we found that for one- dimensional expansion on average large DCC 
domain can not be formed. However, in some particular orbit large instabilities can occur. 
This possibilty also reduces with introduction of friction. However, if the friction is large, 
the system may be overdamped and then there is a possibilty of DCC domain like formation. 
Present paper is an extension of the above study, the spatial part, which had been integrated 
out in our earlier analysis is being studied here. 

Appropriate coordinates for heavy ion collisions are the proper time (r) and the rapidity 
(Y). For a d-dimensional expansion, the change can be effected by the following replacement, 

d d d__ l_d^ 



dt'^ dz^ T dr dr dY"^ 

It can be seen from the above equation that with the introduction of proper time and 
rapidity, a dissipative term comes into effect in the equation of motion. To illustrate the 
role of friction, we further introduce a dissipative term {rj) in the equation. To simplify 
our calculation, we assume boost-invariance as well as cylindirical symmetry in the system. 
With boost-invariance and cylindirical symmetry, eq.^ can be written as, 

— 77-0 — + - + m — -7 = — ^-^ — + 7^ + A w - < > -T /2 < > +Hn^ 

OT^ T OT or'^ r or 

(7) 

where rj is the friction coefficient and we have replaced the fluctuation terms by their coun- 
terpart (T^/2) in the flnite temperature fleld theory |1^,|18[. The equation of motion of 



fields then depend sensitively on the initial temperature of the system. In the following, we 
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assume initial temperature to be Tc = y2/2 — 2ml/X=123 MeV at the initial time rj= 1 
fm. The cooling of the system is described by the equation, 

T d , , 

- + — = 8 

In the weak coupling limit, the friction rj is related to the on-shell plasmon damping 
rate, r] = 2'~fpi. In the standard (f)^ theory, the plasmon damping rate can be calculated for 



coefficient can be obtained as 15 



the a and the 11 fields Assuming that the meson masses are the same, the friction 

V = 27p. = ^s^'-fsAl - (9) 



where fsp = - /f dt^ is the Spence function. At T = = ^2/^ - 2m2/A=123 MeV, 
and if m/T ~ 1, the friction t] = 2.2/m^^. In order to have a comparison, we also consider 
scenarios with ri=0 and ri=0.5 fm~^. We also neglect the temperature dependence of the 
friction coefficient. 

We solve the set of partial differential equations |^ with the quenched initial condition, 
for two dimensional (d=2) expansion using the standard Leap frog method. Evolution of 
the fields upto r=7 fm are followed. As argued before, DCC formation depend critically on 
the initial condition. To refiect the uncertainty in the initial condition, the initial fields are 
randomly distributed to a Gaussian form with the following parameters |T^ , 



<a>=(l-/(r))A (10) 

< TT, >= (11) 
<a^ > - <a >2=< Trf > - < tt, >^= v^/6f{r) (12) 

< a >= < TT, >= (13) 
< >=< TT >2= 4wV6 (14) 

The interpolation function 

/(r) = [1 + exp{r - ro)/T)]-' (15) 

separates the central region where the initial field configuration is different from their vaccum 
expectation value. We use ro=5 fm and r=0.5 fm. 

The equation of motion was solved for 500 trajectories. For each trajectories, at each 
space-time we compute the effective mass ^Iff- The phenomenon of long wavelength DCC 
amplification will occur whenever the effective mass squared is negative. To single out the 
trajectories for which maximum instabilities occur for each trajectories, we calculate the 
following quantity, 

G = y \Tneff\Q{—mlj:j)dTdr (16) 

This can be a measure of instability in a particlar evolution. We call it amplification 
factor. This is an important parameter, as it directly relates to the size of DCC domains. 
In 
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FIG. 1. Distribution of the amplification factor G for the three scenarios (i) /?=0, (ii) r/=0.5 
fm^^ and (iii) r/=2.2 fm~^. 

fig.l, we have shown the distribution of the amphfication factor G for the 500 trajectories. 
Three scenarios; (i) friction free case {ri=0), (ii) small friction (?7=0.5 fm~^) and (iii) large 
friction {ri=2.2 fm~^) are shown. Several intersting features emerges. Maximum amplifi- 
cation factor is increased with friction. However the increase is small {Gmax= 2.19 for ri=0 
and Gmax= 2.55 for r]=2.2 fm~^). Minimum amplification factor consequently decreases 
with friction. As a result, distribution is broadened with friction. The result suggest that 
with friction, while number of orbits with minimum instabilities are increased, some orbits 
can be found with enhanced instabilities. Thus with friction possiblity of DCC domain like 
formation is increased in some orbits. This is contrary to the naive expectation that friction 
will inhibit DCC like formation. However, this is in agreement with our earlier study that 
with friction some trajectories can show appreciable instabilities |jl6[ . 

In fig. 2, we have shown the contour plot of sgn{mf.ff) obtained by averaging the fields over 
the 500 events, as obtained in the three scenarios, (i)?7=0, (ii)?7=0.5 fm"^ and (iii)?7=2.2 
fm~^. Friction lowers the magnitude of '^e// evident from the figure. However, mljrj: 
remain positive throughout the space-time region in all the three scenarios. The result 
indicate that on averaging, quenched initial condition do not lead to DCC domain formation. 
Though not shown, average fields are also found to be order of magnitude lower than obtained 
in an individual event. The result suggests that on average DCC like phenomena will not 
be obserevd. 

The behavior is changed in individual events. In fig.3, we have shown the contour 
plot of sgn{meff) for the most unstable orbit i.e. for which the amplification factor is 
the maximum. The three cases are labeled appropriately. For the friction free case, the 
effective mass remains negative over a considerable space-time. Thus from from r=2-4 fm 
over a considerable radial range (2-4 fm), "^e// is negative and one may expect DCC like 
formation correspondingly. With small friction, the pattern remains essentially the same, 
but space-time region where w^e// remains negative is increased. Entirely different behavior 
results when friction is high {ri=2.2 fm~^). In some radial ranges, soon after the start of the 
evolution, rn^ff becomes negative and remains negative throughout the evolution. This is an 
indiaction of overdamped system. The viscous drag is large and once the system enters the 
unstable region viscous drag forces it remain there. This result is interesting, as it suggest 
possiblity of forming large domain like structure with high viscous force in the system. 
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In fig. 4, the contour plot of tlie tti field for the most unstable orbit are shown. The fields 
have been scaled by For the friction free case several regions can be identified where 
growth of negative tti is considerable. Thus for example around r=6 fm, negative tti grows 
for about 4-5 fm of (proper) time. The results for small friction case is nearly similar. We 
observe that for friction free and small friction case, DCC domain like structure can grow for 
about 4 fm. However, the radial extension is not large, it is about 1 fm. For large friction 
case (?7=2.2 fm~^, the pattern is again different and reminiscent of the nf^^. We find strips 
of positive and negative tti again of radial extension of about 1 fm starts growing just after 
the start of the evolution. Evidently, with large friction the domains can grow for a longer 
duration. 




T fm 



FIG. 2. Contour plot of evolution of sgn{m(,ff) obtained by averaging the fields over the 500 
events for {i)r]=0, (ii)ry=0.5 fm^^ and (m)ry=2.2 fm~^. 
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FIG. 3. Contour plot of evolution of sgn{meff) for the event for which the amplification factor 
is (G) the maximum. Three scenarios (i)?7=0, (ii)?7=0.5 fm~^ and {m)r]=2.2 fm^^ are shown. 

Before we summarise, it is important to note one apparent contradiction we observe 
presently. For the friction free and small friction case domain like structure of tti field do not 
exactly follow the pattern of sgn{meff). Thus in those space-time region where ?^e// is most 
negative, we donot find large amplitude of tti field. Though not shown, we have checked that 
this is true for and tts fields also. However in the large friction scenario, tti follows nearly 
the pattern of "^e//- The reason may be the symmetry breaking term. The assumption 
that phenomenon of long wavelength oscillation occur whenever m^jj is negative implicitly 
neglects the symmtery breaking term. However, for small friction, symmetry breaking term 
may not be negligible. Then, long wavelength oscillations will not follow exactly sgn{meff). 
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FIG. 4. Contour plot of evolution of tti//^ for the event for which the amplification factor (G) 
is the maximum. Three scenarios (i)r/=0, (ii)r/=0.5 fm~^ and {m)r]=2.2 fm~^ are shown. 

To summarise, we have investigated the effect of friction on the possible DCC domain 
formation. In the equation of mation for the linear sigma model fields, we include a friction 
term and solve it assuming boost invariance and cylinderical symmetry. Initial field configu- 
ration was assumed to be Gaussian random with quenched condition (< >= < > = 0). 
At each space-time ^e// calculated. Phenomenon of long wavelength amplification oc- 
curs when it becomes negative, ''^e// obtained by averaging the field over 500 events remain 
positive through out the space-time region investigated, indicating that on average, DCC 
domain like formation is not possible. However, instabilities can occur in some particular 
events. In those events, where the instabilites are the largest, we find evidence of domain 
like structure formation. For friction free and small friction, domains can grow for about 3-4 
fm. But with large friction, the system become overdamped and instabilities can continue 
for longer duration. 
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